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ABSTRACT  

The  design  and  operation  of  a  special  purpose  digital  signal 
processing  system  is  explained.  The  system  consists  of  a  signal 
analyzer,  an  interface  between  the  analyzer  and  the  M.I.T.  Compat- 
ible Time  Sharing  System,  and  appropriate  computer  programs.  With 
inputs  from  an  instrumented  shaft  in  the  M.I.T.  variable-pressure 
propeller  tunnel,  the  system  permits  rapid  spectral  analysis  of 
propeller  generated  forces  and  moments  in  six  degrees  of  freedom. 
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1.  INTRODUCTION 

1.1  BACKGROUND.  On  July  ±>  19635  the  Department  of  Naval  Archi- 
tecture and  Marine  Engineering  at  M.I.T.  submitted  a  prooposal  to  the 
National  Science  Foundation  to  build  a  system  that  would  permit  un- 
steady hydrodynamic  measurements  in  the  variable-pressure  water  tunnel 
operated  by  the  department  (l).  research  into  means  of  reducing  vib- 
ratory excitation  to  the  hull  of  a  ship  depends  upon  the  availability 
of  methods  to  identify  and  measure  periodic  signals  of  known  frequency, 
periodic  signals  of  unknown  frequency,  and  random  signals.  Of  parti- 
cular interest  are  the  six  components  of  periodic  forces  and  moments 
on  the  propeller  shaft  due  to  non-uniform  flow  about  the  propeller  and 
the  fact  that  the  propeller  has  a  finite  number  of  blades. 


The  design  and  construction  of  a  specially  instrumented  shaft 
which  will  enable  simultaneous  measurement  of  the  instantaneous  values 
of  the  three  forces  and  three  moments  due  to  the  propeller  interaction 
is  covered  in  the  Naval  Engineer's  thesis  by  Charles  0.  Norton  (2), 
a  matrix  of  miniature  strain  gauges  and  integrated  circuit  operational 
amplifiers  built  inside  the  shaft  provides  six  signals  in  the  100  mv. 
to  1  volt  range  at  an  impedance  of  a  few  ohms.  All  the  signals  of 
interest  will  be  below  1  khz.  jsaeh  of  the  six  forces  and  moments  to 
be  studied  is  a  linear  combination  of  these  six  signals,  and  the  mat- 
rix of  coefficients  for  this  computation  is  to  be  determined  by  dyna- 
mic calibration  of  the  shaft.  A  separate  circuit  in  the  shaft  housing 
delivers  a  sharp  synchronizing  pulse  once  per  revolution,  jpower  to 


the  amplifiers  and  output  signals  from  the  amplifiers  are  transferred 
between  the  shaft  and  the  outside  via  low  noise  slip  rings. 

The  design  of  a  simple  digital  data-processing  system  to  rapidly 
analyze  the  signals  received  from  the  instrumented  shaft  is  the  sub- 
ject of  this  report, 

1.2  CHuICiJ  OF  SYSTEM.  The  principal  requirement  for  the  data- 
processing  system  were  wide  flexibility  in  types  of  analysis  available 
to  the  experimenter ,  short  turn-around  time,  and  cost  below  $15>,000. 

Originally  it  was  felt  that  all  signal-processing  equipment 
should  be  self-contained  and  located  in  the  propeller  tunnel  experi- 
mental area.  Complete  signal-processing  systems  which  could  correlate, 
signal  average,  and  spectrum  analyze  turned  out  to  be  too  expensive. 
Considered  next  was  a  method  of  recording  data  from  the  experiment, 
converting  it  to  digital  form,  and  then  punching  it  on  paper  tape  or 
cards  for  submittal  to  the  M.I.T.  Information  Processing  Center  for 
"batch"  processing,  -The  long  turn-around  time  made  this  prospect  dis- 
couraging. The  final  choice  was  to  use  a  commercial  digital  signal 
analyzer  in  the  experimental  area,  and  to  make  an  interface  between 
the  analyzer  and  the  IBM  70°U  computer  used  in  the  M.I.T.  .Compatible 
Time  Sharing  System  (CTSS) ,  The  analyzer  would  collect  data  from  the 
experiment,  make  the  analog-to-digital  conversion,  and  either  store 
the  data  from  a  single  sweep  for  analysis  by  an  autocorrelating  or 
crosscorrelating  program  in  the  70°)|,  or  make  hundreds  of  sweeps  in 
averaging  mode,  leaving  in  its  memory  a  digitized  signal  of  enhanced 


signal~to«noise  ratio  for  use  in  a  digital  Fourier  analysis. 

The  signal  analyzer  chosen  was  the  Northern  Scientific  Model  5>5>0, 
with  optional  specifications  as  follows: 

Memory:  102li.  addresses;  10  counts  per  address 

Number  of  inputs:     8  (with  NS-30?  eight  input  expander) 
Readout:  Teletype  Model  33ASKTC  (includes  tape 

punch  and  reader) 

The  NS«5>f>0  rost  the  requirements  of  this  project  in  the  following 
specific  ways: 

1.  Because  it  could  read  its  output  into  a  Teletype,  which  in  turn 
can  be  connected  to  Ci'SS  with  a  simple  electrical  interface  to  a  uata- 
Hione,  the  computer  interfacing  problem  would  not  be  complicated. 

"2*  The  system,  in  the  worst  case,  should  resolve  the  third  harmonic 
of  a  9~bladed  propel] er  operating  at  2000  r.p.m.  with  all  8  inputs 
simultaneously.  With  time  multiplexing  of  the  eight  inputs,  this 
requires  a  sampling  frequency  of  at  least  lU.U  khz,,  or  .-:  a -.sampling 
period  of  no  greater  than  70  microseconds.  The  1JS-550  can  sample  at 
intervals  of  5>0  microseconds. 

3.  The  total  price  of  the  analyzer  and  peripheral  gear  was  about 
$12,000, 

The  electrical  interface  between  the  Teletype- machine  and  the 
Data-Phone  was  designed  and  built  by  the  author.  It  is  electronically 
very  simple  and  requires  no  attention  or  adjustment. 
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The  system  is  completed,  by  a  conversational  program  written  in 
Fortran  IV  which  receives  typed  input  from  the  experimenter  and  the 
signal  analyzer  and  returns  a  Fourier  analysis  of  each  of  the  six 
forces  and  moments.  For  an  analysis  involving  all  1021;  memory  loca- 
tions in  the  analyzer s  the  turn-around  time  is  about  20  minutes s  most 
of  which  is  analyzer  read-out  time.  The  program  listing  and  documen- 
tation is  given  in  Appendix  A. 

.fr'ig.  1  is  a  block  diagram  of  the  completed  signal  analyzing 
system. 


2.  MODIFICATIONS  TO  EQUIPMENT 

2.1  INPUx  MODIFICATIONS.  The  Ntf-550  analyzer  has  an  input  impedance 
of  at  least  100K  and  switch-selected  AC  or  DC  coupling.  The  torque 
and  thrust  outputs  from  the  shaft  will  have  d.c.  levels  considerably 
higher  than  twice  the  a.c.  signal  components,  and  consequently  the 
baseline  adjustments  on  the  analyzer  will  be  insufficient  to  remove 
them.  Unfortunately,  the  time  constant  of  the  a.c.  coupling  is  .01 
second, .  which  is  too  small  to  pass  low-frequency  signals  from  the  ex- 
periment without  significant  attenuation  and  phase  shift.  This  con- 
dition can  be  corrected  by  replacing  the  .1  ]if .  coupling  capacitors  in 
the  eight  input  circuits  with  10  nf .  capacitors  of  equal  voltage  ra- 
tings . 

2.2  SWEEP  MODIFICATIONS.  With  8  inputs  and  a  memory  cycle  time  of 
50  microseconds,  it  takes  ,k  msec,  to  add  one  data  point  to  each  of 
the  eight  memory  subdivisions .  For  simplicity,  assume  that  each  chan- 
nel uses  12^  memory  locations.  We  desire  to  store  data  representing 
one  complete  propeller  shaft  revolution  in  each  subdivision.  The 
maximum  period  of  one  revolution  is  then  123>  X  Ji'msec,  which  corres- 
ponds to  a  minimum  shaft  speed  of  1200  r.p.m.  Increasing  the  speed 
does  nothing  more  than  reduce  the  number  of  data  points  corresponding 
to  one  period. 

Thus,  for  speeds  of  1200  r.p.m.  and  above,  the  time  base  of  the 
NS-550  must  be  set  for  ,k  ms ./point.  Similarly,  for  speeds  between 
600  and  1200  r.p.m.,  a  setting  of  .8  .ms. /point  is  required.  However, 
as  supplied  by  the  manufacturer,  the  NS-5>£0  skips  over  this  range, 
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going  from  .2  ms. /point  to  2  ms. /point. 


The  required  modifications  to  the  tine  base  circuitry  are  not 

difficult,  and  circuit  diagrams  showing  these  changes  can  be  obtained 

2 
from  Northern  Scientific e 

2.3  LINE  FEED  MODIFICATION^ .  Whenever  CT3S  is  called  up  from  a  Tele- 
type terminal,  the  computer  begins  by  signaling  the  automatic  answer- 
back device  on  the  Teletype.  The  codes  which  are  returned  identify 
the  particular  Teletype  machine  and  tell  the  computer  if  it  must  sup- 
ply separate  line  feed  signals  whenever  a  carriage  return  is  either 
sent  or  received.  The  two  codes  are  identical  except  that  in  one  a 
letter  "C"  appears  in  the  two  places  which  contain  the  letter  "D"  in 
the  other  code.  If  the  answer-back  device  is  set  up  for  the  "D"  code, 
the  computer  supplies  a  separate  line  feed.  If  it  is  set  up  for  the 
"C"  code,  the  computer  does  not  supply  the  separate  line  feed,  but  as- 
sumes that  whenever  a  carriage  return  is  either  sent  or  received,  the 
Teletype  machine  will  line  feed  on  its  own. 

The  "D"  code  works  satisfactorily  for  all  input  and  output  ex- 
cept for  signal  analyzer  readout,  because  the  signal  analyzer  produces 
a  line  feed  of  its  own,  and  does  not  wait  for  the  receipt  of  the  com- 
puter-generated line  feed  before  continuing.  This  causes  confusion 
and  loss  of  data. 

Using  the  "C"  code  solves  the  above  problem  and  introduces  ano- 
ther: the  operator  must  type  his  own  carriage  returns  when  typing 
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input  to  the  computer,  and  must  manually  advance  the  paper  Tor  each 
successive  line  of  data  returned  by  the  computer. 

The  solution  to  this  problem  is  to  set  up  the  function  box  in  the 
Teletype  machine  for  "Suppress  Line  Feed"  and  "Hew  Line,  Line  £'eed," 
in  accordance  with  the  technical  manual  supplied  with  the  machine  (3; . 
The  ''Suppress  Line  Feed"  feature  ignores  all  line  feed  signals  re- 
ceived, and  the  "Hew  Line,  Line  .teed"  function  automatically  produ- 
ces a  line  feed  whenever  a  carriage  return  signal  is  either  transmit- 
ted or  received,  uf  course,  the  "C"  answer-back  code  is  used. 
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3.  ELECTRICAL  INTERFACE  BETWEEN  TELETYPE  AlID  DATA-PHONE 

3.1  TELETYPE  SIGNALS,  ihe  33ASRTC  Teletype  that  Northern  Scientific 
supplies  as  an  optional  readout  device  for  the  NS-550  comes  wired  by 
them  for  three  modes  of  operation;  "Tape  to  Type,"  "Analyzer  Head 
•In,"  and  "Analyzer  Read  Out."  In  the  first  mode,  the  tape  reader  sup- 
plies data  to  the  printer.  In  the  second  mode,  the  tape  reader  sup- 
plies data  only  to  the  memory  of  the  analyzer.  The  last  mode  has  the 
printer  (and  tape  punch,  if  desired)  print  (punch)  out  the  data  in  the 
analyzer's  memory.  In  all  these  modes,  the  data  is  transferred  in 
parallel  to  and  from  the  analyzer.  The  parallel-to-serial  conversion 
of  the  code  for  each  symbol  (required  for  transmission  via  the  Data- 
Phone)  is  accomplished  electromechanically  by  the  printer  mechanism. 

The  Teletype  itself  has  provision  for  both  "local"  and  "line" 
operation.  In  local  mode,  an  internal  power  supply  provides  the  60 
ma.  loop  current  required  to  operate  the  selector  magnet  system,  when- 
ever the  machine  is  being  used  for  purposes  not  requiring  a  connection 
to  the  computer,  it  should  be  operated  in  local  mode.  In  line  mode, 
an  external  power  supply  is  required,  and  it  is  connected  to  terminals 
6  and  7  of  terminal  strip  X.  (For  details  of  the  Teletype  wiring,  see 
reference  (k) -  If  the  current  from  the  external  supply  is  modulated 
on  and  off  by  the  ASCII  code  for  a  letter,  the  machine  will  print  the 
letter  (and  punch  it  as  well,  if  desired) .  Alternatively,  if  a  key  on 
the  keyboard  is  -deppressed,  the  Teletype  itself  modulates  the  line 
current  with  the  ASCII  code.   This  also  happens  when  a  symbol  is 
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Fig.  2:  Schematic  diagram  0f  the  electrical  interface 


VJ 

■2. 


•J. 
"2 

w 

-i 

O 

l-^ 

i 

^ 

v> 

V) 

EH 

</> 

£ 

a 

EH 

«0 

CO 

<->N 

+1 

•»    H 

H 

CO     1A 

o 

«» 

-P     XA 

> 

CO 

H     XA 

-p 

O      £5 

■LA 

H 

!>     cm 

CM 

o 

^_^ 

!> 

O 

1A        • 

o 

XA 

H      g 

G 
O 

CM 

O  O 

O 

OvO 

•> 

s — - 

Q)  ^ 

o 

k: 

-3 

O  vQ 

*       o 

& 

CM 
H 

>       CM 
H 

*     H 

.ct 

•»  -^ 

s 

n 

£5 

0-j      s 

o 

CO 

£5      CM 

o 

v— ' 

3-,      -w 

a    o 

r-i 

o 

«    q 

o    v\ 

•H 

r-l 

O      H 

o 

CO 

o 

•H 

•H 

H          <D 

«• 

O 

H          O 

•H         «H 

fO 

^H 

•H         «M 

CO      .3 

O1 

.3 

co    ra 

•  t 

•> 

<• 

H 

CM 

_=t 

GP 

<y 

Z3 

% 

o 

3 


■d 

I 


© 

to 
o 

H 

I 

o 


-!       10 


o 

H 

bO 
bO 
O 
-P 

EH 

Ph 
CO 


u 

CO 


o 

a 

•P  co 
o  S 

•<-!  _ 

-P  xi 

r  !  -rl 
o  fc 
CO  EH 


H     H     H     H 

Q      h     CO      EH 


3 

+i 

u 
o 

*!> 

rt  h 

CO     O, 
to    P, 

^1  to 

°£ 

CO    O 
P-.  f-: 

•P 

fll 

-P    O, 

o    • 

On  O 
I      • 

ft,  xJ 

t)   co 
a  -p 

•H  H 

U    O 
EH    t> 


CM 
EH 


lh 


printed  during  operation  in  "Taps  to  Type"  or  in  "Analyzer  Read  Out" 
mode. 

3.2  DATA-rHONE  SIGNALS.  The  Data-Phone  transmits  and  receives  all 
data  serially >  in  the  form  of  pulses  of  (nominally)  plus  and  minus  1$ 
volts.  Full  details  are  given  in  the  manufacturer's  interface  speci- 
fication (5 J . 

3.3  THE  INTERFACE.  The  interface  translates  the  modulated  line  cur- 
rent from  the  Teletype  to  positive  and  negative  voltages  for  the  Data- 
Phone,  end  vice  versa.  Fig.  2  is  a  schematic  diagram  of  the  unit. 
The  interface  consists  of  a  60  ma.  loop  supply,  a ±10  volt  supply,  and 
two  transistorized  switches.  The  incoming  signal  appears  at  terminal 
3  of  the  "customer's  equipment"  connector.  "Whenever  the  "Mark"  (nega- 
tive) signal  is  received,  the  loop  current  is  switched  on  by  Ql,  -The 
"Space"  (positive)  signal  causes  Ql  to  switch  off.  The  germanium 
diode  connected  to  the  ba.se  of  y2  serves  two  purposes.  First,  its 
back  resistance  is  sufficient  to  assure  that  Q2  is  cut  off  whenever 
the  Data-Phone  is  disconnected.  Second,  it  prevents  damage  to  the 
base-to-emitter  junction  in  the  event  the  Data-Phone  produces  a  nega- 
tive voltage  greater  than  -10  at  terminal  3«  (it  may  go  as  high  as 
-25,  according  to  the  interface  specification.) 

Q3  and  q);  supply  the  outgoing  signal.  The  voltage  divided  con- 
nected to  the  bases  of  Q3  and  Qh  is  chosen  so  that  when  the  Teletype 
is  marking  (terminals  6  and  7  of  strip  X  shorbed)  a  negative  voltage 


1*. 


at  low  impedance  is  presented  at  pin  2}  and  when  the  Teletype  is  spa- 
cing, a  corresponding  positive  voltage  atvlow  impedance  appears ^ at ~ pin 
2. 

The  continuous  positive  voltage  connected  to  pin  20  enables  the 
Data-Phone  to  enter  and  remain  in  the  "Data"  mode. 
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ii.  PROGRAMMING  REQUIREMENTS 

U.l  OUTPUT  FORMAT  FOR  THE  NS-#0.  A  typical  line  of  output  from  the 
NS-5>5>0  consists  of  a  line  feed  (LF),  a  carriage  return  (CR),  and  eight 
6-digit  numbers,  each  separated  by  a  single  space  from  the  preceeding. 
If  each  memory  location  contains  only  5  decades  (which  is  the  case  for 
the  MS-550  used  in  these  experiments),  the  first  digit  of  each  number 
is  always  zero.  The  analyzer  never  prints  a  plus  or  minus  sign.  With 
the  modification  described  in  section  2.3  of  this  report,  the  Teletype 
ignores  the  line  feed,  but  automatically  produces  one  when  it  receives 
the  carriage  return. 

The  initial  CR  and  LF  is  interpreted  by  the  computer  as  a  blank 
card  linage,  and  in  order  to  avoid  having  the  first  8  locations  of  the 
input  array  set  equal  to  zero,  the  program  must  be  written  to  read  the 
blank  image  before  reading  the  first  8  data  points. 

The  following  program  segment  will  successfully  read  JMA.X  points 
from  the  signal  analyzer: 

READ  100 

100  FORMAT  (IX) 

DO  10  NsljJMaX,  8 
NA=N+7 
10     READ  101,  (INDATA(NZ),NZ*N,NA) 

101  FORMAT  (8(I6,1X)) 
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U..2  ALLOCATION  OF  STORAGE  IN  THE  NS-#0.  The  102)4  memory  locations 
in  the  analyzer  are  divided  equally  among  the  input  channels  being 
used.  The  first  end  last  locations  in  each  block  of  memory  are  always 
set  to  zero >  except  for  the  first  location  in  the  first  block,  which 
always  contains  the  number  of  sweeps  made  in  the  experiment.  uata  is 
always  read  out  sequentially  by  memory  address  number.  Data  from  Chan- 
nel A  appears  in  the  first  block  of  memory,  data  from  channel  B  ap- 
pears in  the  second  block,  and  so  on.  A  computed  GO  TO  statement  with 
its  argument  equal  to  the  number  of  input  channels  can  be  used  to  as- 
sign the  data  received  from  the  analyzer  readout  to  arrays  indexed  by 
channel  number  (See  Appendix  A) . 

it. 3  INPUTS  sUr-PLIED  BY  THE  EXPERIMENTER.  Because  time-sharing  allows 
conversational  interaction  between  the  program  and  the  experimenter, 
the  progr?i,i  can  ask  for  the  values  of  various  constants  and  receive 
the  typed  responses  directly.  The  following  sets  of  constants  are  re« 
quired: 

1.  The  number  of  input  channels  used  in  the  experiment 

2.  The  settings  of  the  input  attenuators  on  the  analyzer 

3.  The  matrix  of  coefficients  relating  the  actual  forces  and  moments 

3 

to  the  signals  appearing  at  the  input  of  the  analyzer 

U.  If  the  experiment  involved  signal  averaging,  and  the  readout  did 
not  begin  with  the  first  address  of  memory:   the  number  of  sweeps 
made  by  the  analyzer 

5>.  If  the  program  does  a  Fourier  analysis:  the  number  of  data  points 
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in  a  channel  which  represent  one  complete  shaft  revolution,  and 
the  number  of  harmonics  to  be  computed 

6,  If  the  program  computes  a  correlation  function:  the  range  of  the 
delay  time  (the  argument  of  the  correlation  function) 

7.  If  the  program  crosscorrelates:  the  values  of  the  second  function 
or  ■  an  algorithm  for  computing  these  values 

The  Fourier  analyzing  program  in  Appendix  A  has  numerous  examp3.es 
of  conversational  input-output  routines. 


19 


$.     OPERATING  INSTRUCTIONS  FOR  A  SIX-DIMENSIONAL  ANaLYoIS 

It  is  assumed  in  these  instructions  that  the  reader  is  familiar 
with  Mr.  Norton's  thesis  (2)  which  describes  the  overall  design  of  the 
complete  experimental  system,  the  operation  of  the  NS~5>5>0,  which  is 
covered  in  the  instruction  manual  (6)5  and  the  GTSS  procedures s  which 
are  covered  in  the  CTSS  User's  Manual  (7)* 

PROCEDURE 

5.1  Push  the  STOP  button  on  the  NS-550. 

5.2  Set:  DECADE  SELECT  to  1CT 

READ  MODE  to  CRT 
aDD~SUB  to  ADD 

MEMORY  CROUP  to  1/1  - 

TIME  BASS  to  2  ms. /point 

NO.  OF  SWEEPS  T0<*> 

DELAY  TIME  to  1  ms. 

PRE  DELAY/  POST  DELAY  to  OFF 

ADDRESS  SELECT  to  OFF 

VERTICAL  OFFSET  to  OFF  • 

CONVERSION  GAIN  tu  kVE 

SWEEP  MODE  to  AVE 

TRIGGER  to  RECURRENT 

VIEW  INPUT-VIEW  MEMORY  to  VIEW  MEMuRY 

INPUT  SELECT  to  EXT 
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5.3  Push  START  READOUT,  k  pattern  of  some  type  should  appear  on  the 
CRT.  xf  nothing  appears,  carefully  adjust  the  INTENSITY,  HORI- 
ZONTAL, and  VERTICAL  controls  until  something  does.  Always  keep 
the  INTENSITY  control  set  to  the  minimum  level  which  produces  a 
visible  trace  or  spot  on  the  CRT,  This  is  especially  important 
when  only  a  single  spot  appears  on  the  CRT. 

$M  Push  the  ERAS.&  button.  The  CRT  pattern  should  immediately  degen- 
erate to  a  straight  line.  Adjust  the  VERTICAL  POSITION  control 
so  this  line  appears  at  the  very  bottom  edge  of  the  CRT  face. 
Adjust  the  HORIZONTAL  POSITION  and  HORIZuNTAL  EXPAND  controls  so 
that  this  line  extends  slightly  less  than  the  width  of  the  CRT 
face. 

£.£  Push  STOP.  Then  push  START  MEASURE.  Let  the  unit  run  for  about 
3-0  seconds,  then  push  STOP,  followed  by  START  READOUT.  This  cas- 
ses  a  small  amount  of  noise  to  be  read  into  the  memory. 

5.6  Set  the  DECADE  SELECT  to  10  and  observe  that  10  horizontal  lines 
appear  on  the  CRT.  Using  the  VERTICAL  EXPAND  and  VERTICAL  POSI- 
TION controls,  make  these  10  lines  coincide  with  10  of  the  11 
scribed  reticle  lines.  This  completes  calibration  of  the  CRT. 

5.7  Push  STOP.  Set  DECADE  SELECT  to  10  .  Set  the  INPUT  ATTENUATuRS 
for  all  k  inputs  of  the  NS-550  "and  all  i>  inputs  of  the  NS-305 
to  1  volt.  Set  all  8  of  the  AC-DC  switches  to  AC. 

5.8  Connect  the  six  strain-gauge  outputs  from  the  shaft  to  channels 
A  through  F.  Connect  the  synchronizing  pulse  cable  from  the 
shaft  to  the  positive  input  of  Channel  H  and  also  to  the  EXT 


21 


TRIGGER  connector* 

5.9  Start  up  the  propeller  and  set  the  shaft  r.p.m.  to  the  value  that 
will  be  used  in  the  experiment.  Set  the  THE  BASE  to  »k  ms./ 
point  if  the  r.p.ra.  is  above  1200.  Set  the  TIME  BASE  to  .8  ms./ 
point  if  the  speed  is  between  600  and  1200  r.p.m.  Set  the  TIME 
BASE  of  the  NS-305  to  the  same  value  that  is  set  on  the  NS-550, 

5.10  Set  VIEW  INPUT  to  A.  Set  TRIGGER  to  -H3XT,  and  push  START  MEASURE. 
Adjust  the  TRIGGER  SENSITIVITY  for  stable  triggering  of  the  hori- 
zontal sweep,  as  indicated  by  a  steady  pattern  on  the  CRT  that 
doss  not  move  left  or  right.  Center  the  pattern  vertically  by 
using  the  Channel  A  BASELINE  ADJUST.  If  the  pattern  goes  off  the 
top  or  bottom  of  the  CRT,  set  the  Channel  A  INPUT  ATTENUATOR  to 
10  volts.  If  the  pattern  is  very  small,  set  the  Channel  A  INPUT 
ATTENUATuR  to  0.1  volt. 

5.11  Switch  VIEW  INPUT  to  u.  Using  the  Channel  B  BASELINE  ADJUST  and 
INPUT  ATTENUATOR,  position  and  vertically  scale  the  pattern  as 
before,  itepeat  for  Channels  C  and  D. 

£♦12  Set  the  DECIDE  SELECT  to  LOG,  and  VIEW  INPUT-VIEW  MEMORY  to  VIEW 
MEMORY  (This  permits  viewing  of  the  four  inputs  of  the  NS-305) • 
Set  VIEW  INPUT  on  the  N5-305  to  E.  Using  the  Channel  E  INPUT  AT- 
TENUATOR and  BASELINE  ADJUST,  position  and  scale  the  pattern  as 
before.  Repeat  for  Channel  F. 

5.13  set  the  VIEW  IN±UT  to  G.  Center  the  line  using  the  channel  G 
BASELINE  ADJUST. 

5.1U  Set  VIE1.'/  INPUT  to  II .  The  pattern  should  be  a  straight  horizontal 
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line  except  for  a  pulse  near  the  right-hand  side.  With  the  EASE- 
LINE  adjust  for  Channel  H,  position  the  straight  line  about  l£^> 
up  from  the  bottom.  With  the  INPUj.1  ATTENUATOR,  make  the  pulse 
height  at  least  two  divisions  on  the  reticle. 

£.l£  Push  STOP.  Set  DECADE  SELECT  to  10  .  Push  START  READOUT.  Push 
ERASE,  rush  STOP.  This  completes  the  adjustments. 

£.16  To  begin  the  measurement,  push  START  MEASURE.  Initially  the  CRT 
will  show  only  a  horizontal  line  at  zero.  As  the  experiment  con- 
tinues, vertical  deflections  will  appear  and  grow.  Set  DECADE 

£ 

SELECT  to  10  .  Continue  the  experiment  until  the  position  of  any 

portion  of  the  trace  approaches  the  top  of  the  CRT.  Then  push 
STOP.  (Alternative:  Set  the  NO,  OF  SWEEPS  control  to  800  before 
pushing  START  MEASURE.  The  analyzer  will  make  800  sweeps  and 
stop  automatically.)  This  completes  the  measurement  phase  of  the 
experiment . 

£.17  Turn  on  the  electrical  interface  between  the  Teletype  and  the 
Data -Phone,  Turn  on  the  Teletype  in  "Line1'  mode.  Set  the  func- 
tion switch  to  "ANALYZER  READ  OUT."  Dial  up  the  computer  and 
LOGIN.  Then  type  LOADGO  PROP.  The  program  will  begin  by  asking 
for  a  number  of  constants.  The  following  steps  give  the  details. 

£.18  "TYPE  NUMBER  OF  INrUT  CHANNELS."  Format  is  II.  Type  in  the  num- 
ber six.  (The  sync  input  on  Channel  H  is  not.  counted  as  an  in- 
put.) 

£.19  "TYPE  INPUT  ATTENUaTOxt  SETTINGS."  Format  is  6F7.2  Type  in  the 
settings  of  the  six  INPUT  ATTENUATORS,  A  through  F. 
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5.20  "TYPE  NUMBER  0?  OUTPUT  CHANNELS  DESIRtffi."  Format  is  II.  Type  6. 

5.21  "TYPE  NUMBaR  OF  HARMONICS  DESIRED."  Format  is  12.  The  basic 
period  of  analysis  is  one  shaft  revolution,  EXAMPLE:  To  obtain 
analysis  up  through  the  third  harjaonic  of  the  blade  rate  for  a  5- 
bladed  propeller,  1$  harmonics  are  required.  The  maximum  number 
that  can  be  typed  in  is  5>0« 

5.22  "ENTER  THE  N-OUT  BY  N-IN  CALIBRATION  MATRIX."  The  response  must 
be  the  6x6  strain-gauge  calibration  matrix  C,  . ,  where  k,  the 
row  subscript,  is  the  output  channel  number,  and  i,  the  column 
subscript,  is  the  input  channel  number.  The  dimensions  of  each 
entry  must  be  either  ft.-lbf ,/volt  or  lbf ,/volt,  and  the  computer 
will  read  six  successive  rows  of  six  numbers  each  on  a  6i»'7.2  for- 
mat. 

5.23  "TYPE  NUMBER  OF  DATA  POINTS  TO  BE  RR4.D."  Format  is  ik.  Type  in 
the  number  1021;. 

5.2lj  "HERE  IS  HOW  I  HAVE  YOUR  SETTINGS-"  The  computer  will  read  out 
all  the  values  as  it  has  them,  and  gives  self-explanatory  in- 
structions on  what  to  do  if  there  is  an  error  anywhere. 

NOTE:  In  the  nexb  two  instructions,  the  program  mentions  switch- 
ing to  either  manual  or  automatic  line  feed.  If  the  modifica- 
tions described  in  section  2.3  of  this  report  have  been  made  to 
the  Teletype,  the  line  feed  instructions  should  be  deleted  from 
the  program. 
5.25  "I  AM  READY  TO  ACCEPT  DATA  FROM  THE  SIGNAL  AVERAGER."  If  a  paper 
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tape  of  the  readout  is  desired,  turn  on  the  tape  punch  and  hit 
the  space  bar  on  the  keyboard  about  10  times  to  make  a  leader  on 
the  tape,  but  do  not  type  a  carriage  return  or  a  line  feed.  Push 
STOP  on  the  NS-5&0.  Set  the  READ  MODE  to  TYPE  and  push  START 
READOUT.  The  readout  will  take  about  If?  minutes* 

£.26  "SWITCH  TO  AUTOMATIC  LINE  FEED  AND  TYPE  CARRIAGE  RETURN."  See 
the  note  above.  If  this  instruction  and  the  read  statement  im- 
mediately following  it  in  the  program  (Appendix  A)  have  not  been 
deleted,  simply  hit  the  carriage  return  to  continue. 

£.27  "TYPE  NUMBER  OF  DaTA  POINTS  COMPRISING  ONE  CuMPLErE  PERIOD."  In 
the  event  the  program  was  not  able  to  detect  the  sync  pulse  in 
Channel  8,  it  will  make  the  above  request.  If  this  happens,  look 
at  the  last  16  rows  of  the  printed  readout.  Observe  that  the 
sixteenth  row  from  the  bottom  begins  with  000000.  Number  "these 
rows  1  through  16  starting  with  the  one  beginning  with  zeroes. 
Then  find  the  discontinuity  in  the  printed  data  values.  If  R  is 
the  row  and  C  is  the  column  where  the  discontinuity  appears,  the 
number  to  type  is  8(R-1J4-(G«1),  The  format  is  lU« 

£.28  "REaD  OUT  NUMBER  OF  BOSSES."  If  the  computer  receives  a  zero  for 
the  first  number  read  in  from  the  analyzer,  it  trill  make  the 
above  request.  Push  STOP  on  the  Ns-££0.  Set  MEMORY.  GROUP  to  1/1 
.  and  turn  the  ADDRESS  SELECT  switch  OFjf.  Set  READ  MODE  to  -CRT, 
Push  START  READOUT  and  then  push  STuP.  Set  READ  MODE  to  TYPE  and 
push  START  READOUT.  After  the  first  number  is  printed  and  before 
the  eighth  is  printed,  push  STOP,  set  READ  MODE  to  CRT,  and  then 
hit  the  carriage  return. 
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5>.2°  The  computer  i-rill  now  read  out  by  channel  number  and  harmonic 
number  the  amplitude  in  lbf e  or  ft ♦-lbf  and  the  phase  in  deg- 
rees of  each  harmonic.  The  2eroth  harmonic  is  the  average }  and 
it  has  no  significance  in  this  type  of  experiment^  since  tlie  in- 
puts are  all  a.c.  coupled.  The  channel  numbers  can  be  related  to 
a  specific  force  or  moment  by  a  brief  consideration  of  the  compu« 
tation  implied  by  the  6  A  6   calibration  matrix. 

5.30  "THIS  COMPLETES  THE  ANALYSIS.  TYPE  1  TO  START  AuAIN,  0  TO  STOP." 
If  a  "1"  is  typed,  the  program  will  return  to  step  $»?3. 
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I 
APPENDIX  A:   "PROP 5"  A  COMPUTER  PROGRAM  FOR  FOUitlER  ANALYSIS 

The  program  described  in  this  section  embodies  the  principles 
outlined  in  section  h  and  is  the  progrrjn  used  for  the  6-dim3nsional 
analysis  of  section  £.  The  listing  begins  on  page  32.  PROP  can  be 
divided  into  5>  parts: 

1.  TYPED  INPUTS:  PROP  requests  values  from  the  experimenter  and  pro- 
vides an  opportunity  to  start  again  if  errors  have  occurred.  (Cards 
1-37) 

2.  ANaLXZER  READOUT:  The  program  reads  the  analyzer  output  into  the 
array  INDATA,  and  then  assigns  the  values  to  appropriate  blocks  in  the 
array  INPUT,  according  to  the  rule  stated  in  section  lj.,2.  If  5>  or 
more  input  channels  are  being  used,  the  program  scans  the  contents  of 
Channel  8,  looking  for  a  discontinuity  of  at  least  20,000  counts.  If 
it  finds  such  a  discontinuity,  it  sets  the  number  of  data  points  rep- 
resenting one  complete  shaft  revolution  equal  to  the  location  of  the 
discontinuity.  If  the  discontinuity  is  not  found,  or  if  the  number  of 
inputs  is  k  or  less,  PROP  requests  the  experimenter  to  type  in  the 
correct  value.  (Cards  38-8?) 

3.  MATRIX  MULTIPLICATION:  The  program  converts  strain-gauge  signals 
into  forces  and  moments c  (Cards  88-98) 

ii.  FOURIER  ANALYSIS:  PROP  continues  by  calculating  the  coefficients 
for  the  Fourier  sine  and  cosine  series  for  each  force  and  moment. 
(Cards  99-116) 
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5>.  OUTHJT:  The  program  determines  the  number  of  signal  analyzer  pas- 
ses (sweeps),  calculates  the  magnitude  and  phase  of  each  harmonic  for 
each  output  channel,  prints  the  results,  and  provides  the  experimenter 
the  opportunity  to  read  in  new  data  for  the  nexb  analysis.  (Cards  117- 
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The  following  is  a  chart  of  the  variables  used  in  tne  program. 


VARIABLE  HAMS    RANGE  OF  VALUES 


DESCRIPTION 


A(K,L) 


AMP 

ave'; 


B(K,L) 


CALIB(K,I) 
FuLSCA(l) 


biax 


lbf./100  counts  or 
ft.-lbf ./100  counts 

ft.-lbf.  or  Ibf. 
ft.-lbf.  or  Ibf. 

Ibf ,/100  counts  or 
ft.-lbf./lOO  counts 

27T/Ri\|DPC 

Ibf ./volt  or 
ft.-lbf . /volt 

.02^-100.  volts  per 
100  counts 

1-8 

1-8 


Fourier  sine  coefficient  for 
the  Lth  harmonic  in  the  Kth 
output  channel 

The  magnitude  of  a  harmonic 

The  average  value  of  the 
signal  in  an  output  channel 

Fourier  cosine  coefficient 
for  the  Lth  harmonic  in  the 
Kth  output  channel 

Coefficient  used  in  Fourier 
series  calculations 

The  k,i  entry  in  the  strain- 
gauge  calibration  matrix 

Input  attenuator  setting  of 
input  channel  i 

Index  giving  input  channel 
number 

Number  of  input  channels  be- 
ing used  in  the  experiment 
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INDATA(N) 


IwPUT(l,J) 


0-99999 


0-99999 


J,  JP,  dPO,  etc.  1-1022 


JMAX 

K 

L 

M 


NA 

NDPC 


NOUTCH 


PHAS4 


1-1021; 

1-8 

l-£0 

i-5o 

8,  ±6S  2i\.}  etc, 
1-1022 


1-8 


ta9  m 

0-1 

NZ 

1-1021+ 

OUTPUT (K, J) 

lb./lOO  co imts  or 

ft.-lbf.AOO  counts 

PASSES 

1.-99999. 

-90  o  to-V2?0.  degrees 


The  Nth  data  point  read  in 
from  the  signal  analyzer 

The  Jth  data  point  in  input 
channel  i 

Index  giving  data  point  num- 
ber 

Number  of  data  points  read 

Output  channel  index  number 

Index  of  harmonic  bsing  cal- 
culated 

Number  of  harmonics  the  pro- 
gram is  to  calculate  for 
each  output  channel 

Index  used  in  reading  INDATA 

Number  of  data  points  per 
channel  which  represent  one 
complete  period  of  the  input 
signal 

Number  of  output  channels 

Branching  variables 

Index  used  in  reading  INDATA 

The  Jth  data  point  in  output 
channel  k 

The  number  of  sweeps  made  by 
the  signal  analyzer  in  an 
experiment 

The  phase  of  a  harmonic  in 
reference  to  a  sine  function 
of  the  same  frequency 
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RBI 
RJ 

scos(jf~;:~* 

SSIN(J)  *** 


0. -99999 o 
1.-1022. 
-1.0  to  1.0 
-1.0  to  1.0 


Real  value  of  INPUT (I, J) 
Real  value  of  J 
Stored  value  of  a  cosine 
Stored  value  of  a  sine 


"AVE  is  printed  as  the  0th  harmonic  in  the  output  format, 


*-* 


PHASE  is  initially  calculated  in  radians  from  the  arctangent  lib- 
rary routine 3  and  then  is  converted  to  degrees }  taking  account  of 
the  quadrant. 

"SSIN  and  SCOS  are  used  to  avoid  repeated  computations  of  the  sine 
and  cosine  values  needed  for  the  Courier  analysis  of  a  given  har- 
monic in  different  output  channels. 
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APPENDIX  B:  NOTES 

1.  Representing  the  low-impedance  source  as  an  ideal  voltage  source, 
the  transfer  function  for  the  capacitor  coupled  input  is 

urn.  ^Rc* 

A  straightforward  calculation  shows  that  the  phase  of  H(f )  is 
and  the  magnitude  is 


With  the  shaft  r.p.m.  set  to  °60,  f  is  16  hz,  and  with  RC  ,01  sec,, 

Increasing  the  value  of  RC  to  1  second  permits  reducing  the  r.p.m,  to 
only  96,  Under  these  conditions ,  the  phase  shift  is  only  5>,7  degrees 
and  the  magnitude  is  down  by  a  factor  of  only  .005. 

2.  Northern  Scientific,  Inc.,  2551  West  BeltLLne,  P.O.  Box  66,  Mid- 
dleton,  Wisconsin,  53562.  Phone  (608;  836-6511. 

3.  In  the  computer  program  of  Appendix  A,  each  entry  in  the  calibra- 
tion matrix  is  assumed  to  be  independent  of  frequency.  If  dynamic 
calibration  of  the  shaft  shows  that  this  assumption  is  not  valid,  the 
program  will  have  to  be  rewritten  to  do  a  Fourier  analysis  of  the  in- 
coming data  before  doing  the  conversion  of  strain-gauge  signals  to 
force  and  moment  outputs,  xhen  each  harmonic  will  have  to  be  convert- 
ed separately  in  the  program,  using  values  for  the  calibration  matrix 

entries  that  are  appropriate  for  the  frequency  of  the  harmonic. 
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